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1. INTRODUCTION 

The DTC method is based on a direct control of the electromagnetic torque and stator flux of 
the induction machine. This technique was firstly introduced by I. Takahashi and T. Noguchi in the mid 
eighties [1], [2] in order to give an alternative control to the existing vector control. The DTC has excellent 
torque response, but it causes many oscillations due to the control of the flux and the torque through 
hysteresis band controllers, which influence its performance and implementation [3]. For economical and 
technical reasons it is preferable to remove the sensor of the rotor speed. Several strategies have been 
proposed in order to estimate the rotor speed in the IM sensorless drive. Among these techniques, the model 
reference adaptive system (MRAS) is one of the most commonly used techniques mainly to its relative 
simplicity and low computational effort [2], [4], [5], [6]. 

The reactive power MRAS-based is robust against the variation of the stator resistance while 
avoiding pure integration, but suffers from a problem of instability, therefore the rotor flux MRAS-based 
proposed by Schauder [7] is the most popular strategy of MRAS techniques, and much work has been done 
on improving the performance of this strategy. In general, the main problems of the sensorless control are 
operating at low speed, the problems of sensitivity and the variation of parameters of the induction machine, 
specifically, the stator resistance and the rotor-time constant, with also the integration problems [8], [9]. 
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For any estimation-based model, the problems are localized at low speeds, because the voltages induced in 
the rotor are very small and disappear at zero stator frequency [9], and stator voltage drop becomes influent, 
thus the determination of the stator resistance is required for good implementation of the drive. On the other 
hand, the rotor-time constant is required in the MRAS for estimating the rotor speed. These parameters have 
the most significant effect on the accuracy of the speed estimation; their precise knowledge is of paramount 
importance to obtain a good speed control. 

The proportional integral (PI) controllers are widely used in the application of industrial control 
systems. They have a simple structure and can provide a satisfactory performance in a wide operation range. 
Therefore, the most adaptive systems described for observers MRAS acquirer linear PI regulator for 
generating the estimated speed, but due to the variation of the parameters of the induction machine mainly 
related to the heating of the machine and integration drift problems, we have introduced new types of 
controllers based on sliding mode in the speed regulation and artificial neuron network thechnique for 
the conjoint estimation of the most time-varing parameters for his simplisité and the ease of implementation 
compared to the other methods. 

From these methods, the sliding mode control initially investigated for electric motors by Utkin [10] 
has the finest results that lead to the step up of IM control. In the other hand, this type of controllers has a big 
chattering problem [6]. Therfore, a new theory called Reaching Law Control method (RLC) has been 
introduced in this sliding mode controller and many researchers have been conducted in this way. In RLC 
method the error dynamics is specified in the reaching mode, with selecting appropriately the parameters 
both, the dynamic quality of the SMC system can be controlled and the chattering is uninvolved [11]. 

In this paper, we propose to replace the conventional PI controllers used in the speed regulation by 
a new exponential reaching law (ERL) technique based on sliding mode controller. Then, we propose 
a conjoint on-line estimation of the stator resistance and the inverse of rotor time constant by one block of 
artificial neural network technique based on sensorless MRAS-DTC control. 


2. DIFFERENTIAL EQUATIONS OF THE STATE MODEL 

The DTC is built in a steady state. It does not need the stator angle, so the state model equations will 
be written taking into consideration the stator currents Zso, Isg and, the rotor flux a, @-pas state variables as 
shown in (1) 
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Where R, and R, are the stator and rotor resistances, L, and L, are the stator and rotor inductances, Lm 1s the 
. M? 
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LsLr 





. The mechanical equation (2) is: 
dn 
Ja = Te — Tr = fo (2) 


Where J is the moment of inertia, Te and T, are the electromagnetic and load torques, f, 1s the viscous friction 
coefficient, and 2 is the mechanical rotor speed. 


3. THE DTC WITH SLIDING MODE SPEED CONTROL 
3.1. Theory of DTC and sliding mode with ERL 

The direct control of the torque and flux is based on the error between the reference values and 
the estimated values of the torque and flux; it is possible to control the states of the inverter to reduce errors 
directly in the band limits of the predetermined hysteresis controller. The stator flux and the electromagnetic 
torque of the machine are obtained from (3-4): 
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d 
a Ps = Vs — Rs h 
t 
Ps = Pso + Ja Vs — Rs * Is) dt 


Te = Pp (Psalsp i Osplsa (4) 


With p, is the poles pair number. This is to keep the two quantities within the hysteresis band. In addition, 
the voltage vector applied in each switching time is obtained at the output of this regulator with the help of 
the Takahashi Table 1. 

In sliding mode control the trajectory of the system is conducted toward a sliding surface and kept 
on it by means of a suitable switching logic until it reach the equilibrium [10]. The conditions of convergence 
enable the dynamic system to converge to the sliding surfaces S(x). Two conditions are noted wich 
correspond to the reaching mode of the system state [11] as shown in (5). 


S(x)S(x) < 0 (5) 


This is to make a positive scalar function V(x)> 0 for the state variables of the system. This function is 
typically used to ensure the stability of nonlinear systems. Defining the Lyapunov’s function by (6): 


V(x) = 55x) (6) 


To check the reaching condition on the sliding surface, the classical reaching law proposed in sliding mode 
control is presented as shown in (7): 


S = —k sign(S), k >0 (7) 


The time required for the system to attain the sliding surface, called reaching time t, 1s (8): 


[Sol 
t, = e (8) 


As in [12], the suggested reaching law is based on the choice of an exponential expression that adapts with 
the changes of the switching function. This exponential reaching law is presented as shown in (9-10): 


: k . 
S= -yg Seg), k>0 (9) 
N(S) = 0) + (1 — Op )e7 #5 (10) 


0o Is a strongly positive offset smaller than one, p and a are also strictly positive constant. We note that 
the reaching law presented in (9) does not influence the stability of the control process, because N/(S) is 
strictly positive. From (9), we can observe that if |S| increases for a few reasons, N(S) approaches 0o, hence 
k/N(S) converges to k / 05, which is larger than k, this means that k/N(S) increases, thus enhances the pull of 
the sliding surface at the same time. In addition, if |S| decreases, N(S) approach one, so k/N(S) converges to 
k. This can prove that when the system approaches the sliding surface, k/N(S) progressively decreases which 
limits the chattering when the system attains the sliding surface. Therefore, the suggested reaching law 
permits the sliding mode controller to dynamically adapt to the change of the switching function by 
permitting k/N(S ) to vary from k to k/@) . Note that if 05 is selected to be identical to one, the reaching law 
presented in (9) becomes the same reaching law in (7). 


3.2. ERL sliding mode speed controller 


The SM-ERL speed controller presented in Figure 1 which is transient dynamic response depends 
on the selection of the sliding surfaces as shown in (11-12). 


S(t) = ealt) — Jy (ky — a) ea (t)d(t) (11) 
ealt) = A(t) - 2*(t) (12) 
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Where: k; is a negative constant and Q* is the reference rotor speed. 
By derivating (12) and using (2), where f;=T,/J, a=f,/J and b=1 /J one can write (13): 


éo (t) = A(t) — N*(t) = —aeg (t) + u(t) + d(t) (13) 


Where d (t) is the total uncertainties of the parameters a, b, and fı, and the following terms are lumped in u (t) 
as shown in (14): 


u(t) = bT.(t) — aR (t) — f(t) — 2* (0) (14) 


The switching control is chosen as shown in (15-16): 





u(t) = kieo (t) — Ky (S) sign(S) (15) 
Ky(S) = 7% (16) 


With the assumptions that kz defines the upper bound of the uncertainties (k => |d(t)|, XV t), and (k; - a) <0, 
so the speed tracking error eo(t) tends to zero exponentially as the time tends to infinity. The torque current’s 
reference is given by replacing (15) in (14): 


L= b~t|kieg — K,(S)sgn(S) + aQ* + 0* + fil (17) 


Additionally, the problem of chattering is handled by using a soft function (18) in the sliding surface 
neighbourhood, as it’s usually done in the variable structure control: 


=S, IS| <€ 
sign(S), |S| 2 € 


sat(S) = | (18) 


Where the constant factor e defines the thickness of the boundary layer and sat is the saturation function. 
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Figure 1. The Simulink diagram of SM-ERL controller 


4. CONJOINT SPEED AND ONLINE ESTIMATION OF PARAMETERS 
By using the rotor-flux based MRAS, the rotor speed, the stator resistance and the rotor-time 
constant are reconstructed from the reference and the adjustable model [7] presented in Figure 2. 
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Figure 2. Simulink of MRAS thechnique used 


The reference model or the voltage model given by (19, 20) generates the rotor flux. The same model is seen 
as an adjustable model for the online estimation process of the stator resistance [7] [13]. 


Ly x 
PPra = M (kem~ Rdam OL Plu (19) 


Ly ~ 
PPrg = M (Vs, - Rslsg = OL plz) (20) 


In the adaptive model, or the current model, the equation describes the components of the rotor flux 
expressed in terms of stator current components, rotor speed and rotor-time constant. The rotor flux 
components obtained from the adaptive model are given by (21-22) [7]: 


ia M 1 a AJA 

PPra = 7, lsa -— T Pra — Wr Prge (21) 
Pi M 1 a a Aa 

PPrg = 7, ‘sp a 7, Pre + Wr Pra (22) 


The MRAS structure proposed in this paper is based on the estimation of the stator resistance, 
rotor-time constant and the rotor speed which is designed upon the concept of hyper stability [14]. In order to 
make the system asymptotically stable, the proportional-integral controller is used as shown in (23-24): 


Ki 
Wr = Bark + m (23) 


Cas ~ PraPrp > Orp Pra (24) 


The equations of the errors used in the estimation of the stator resistance and the rotor-time constant 
are given by (25-26) [13] [14]: 


Crs Z isa (Pra — Pra) + Isg (Prg — Prp) (25) 
eu = Isa (Gra — Pra) ag Isg (Prp n Prg) (26) 


The adaptation mechanisms which minimize the error between the two models are designed with 
artificial neural networks. In a neural controller type there are three stages, the first is usually the calculation 
of the weighted sum of the inputs [15] [16]. The obtained value is compared to a reference value, which is 
the bias (b). From this value, transfers function calculating the value of the neuron state and the equation as 
shown in (27). This value will be passed downstream neurons where the output of a neuron is always a scalar 
[17] [18]. 
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The network is designed in Matlab (NNtool.) by using a single neural block with three inputs which 
are xi = ls x= (9, -G,,) +(,,—G,,) and x3= Isg for the estimation of two parameters &, R,). 


The back-propagation of the error algorithm is used. It consists of two layers: a hidden layer with five 
neurons, with a hyperbolic tangent activation function (tansig), and an output layer with four neurons, 
the activation function is linear (purelin). The learning method was made initially by gradient descent 
(traingd) and then the Levenberg-Marquardt (trainlm). 

Finally, we conclude the scheme used in the experimental tests using Matlab Simulink and the RTI 
blocks of the SPACE (DS 1104). In Figure 3 is control schema used in simulink and experimental tests. 














MRAS 


Figure 3. Control schema used in simulink and experimental tests 


5. EXPERIMENTAL IMPLEMENTATION 

The experimental test bench presented on the Figure 6 consists of a three-phase Y connected squirrel 
cage induction motor: 1 kW, 2880 rpm, 230/400V, 4/2,8A. The electrical and mechanical parameters are 
identified with the conventional methods which are cited in the Table 1. The motor is loaded by a DC 1kW 
generator with resistive load; the control system is implemented on a dSPACE DS1104 card using 
the Matlab/Simulink RTI block. Sensors with Hall Effect are used to measure the stator phase currents, and 
the DC voltage, the motor is coupled with a 1024 pt incremental encoder for the rotor speed measuring. At 
first, a trapezoidal speed reference is used, changing between -2000 rpm and 2000 rpm, and loaded up to 
50% of rating load by the DC generator. The low-pass filter is a good solution for reducing the chatter of 
the reference torque obtained from the SM controller, but the cut off frequency of this low-pass filter affects 
the performance of the regulators. Using small values reduces torque ripples reference but introduces more 
delay in the system dynamics. A frequency of 60 Hz was found to be a good compromise between torque 
ripples and dynamic response. In Figure 4 is snapshot of the experimental test-bench. The parameters of 
the machine are identified by conventional methods which are cited as following. 


Table 1. The machine parameters 


Variable Value Variable Value 

Number of pole I Rotor inductance 0.7490 H 

Stator resistance 6.58Q Mutual inductance 0.7209 H 

Rotor resistance 5.81Q moment of inertia 0.00207kg.m? 

Stator inductance 0.7490 H viscous friction coefficient 0.000173N.m/(rad/s) 
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Figure 4. Snapshot of the experimental test-bench 


The Figures 5 are organized as follows: from the top down and from left to right are the wave forms 
of the reference (N°) and actual (N) and observed (Noss) rotor speed, the electromagnetic torque (T) the stator 
phase current (isa) and the rotor flux (ọ,). On Figures 6 and 7 the estimated stator resistance (R.), and 


the estimated inverse rotor-time constant ( 1/T,-) are added. 


5.1. Study of the performance of the SM-DTC with MRAS-PI based speed estimation: 

The values of the PI regulator for speed estimator MRAS are: Kp=20, 15 and Ki=7076, 
the sampling times used is Ts =100ms and solved by the Euler method. In Figure 5 is Graphic of the SM- 
DTC with speed estimation MRAS-PI. 

At 50% load (Figure 5.), comparing with [6] where the author used a sliding mode controller for 
the regulation of the estimated speed by a MRAS technique, who found a problem of a small overshoot 
during the change of speed. Therefore, the ERL method integrate the controller sliding mode in this paper 
managed to solved the problem of the speed overshoot with a very good dynamic reply as shown in (5.a) 
where the actual speed (N) exhibits the same behavior of that observed with negligible estimation error 
are follows the trapezoidal reference (up to + 2000 rpm) instantly without overshoot or static error. 
But, a small divergence is seen when passing through zero, which is immediately corrected. The drive is 
stable even in the transient phase with little chattering in the torque (5.b) and the flux (5.c), thanks to the 
ERL-SM controller. 





Nírpm) 


Electromagnetic torque { N.m} 











current,a(A) 





time(s) 


(d) 
Figure 5. Graphic of the SM-DTC with speed estimation MRAS-PI 
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5.2. Study of the performance of SM-DTC with MRAS-PI based speed estimation and with ANN 
parameters estimation 

In Figure 6 the machine is operating with 50% of the rated load. Compared to Figure 5 we see 
the benefit of the online conjoint parameters (R,, 1/T, ) estimation by an artificial neural networks 
controller. First, to test the performences of this control at high speed, we used a trapezoidal reference speed 
varies between 2500 and -2500 rpm. As shown in Figure 6(a) the both estimated and measured speed 
following his reference instantly without any speed drop but with negligible static error that is not affected 
the stability or the dynamic response of the drive. The two Figure 6(b), 6(d) shows compared to 5(b), 5(d) 
a very good improvement in the waveform of the stator flux and electromagnetic torque. It is clear that 
the chattring in these waveforms is remarkably reduced specially during the speed changing, this 
improvement mainly attributable to the proper definition of the estimated parameters by MRAS based ANN 
thechnique, which is very stable, and is not deteriorated even at the zero crossing as shown in Figure 6(e) and 
6(f). Therefore, we can say that our method is one of the most effective Technics in the field of online 
parameters identification and the easier to implement compared to the other methods as in [19] where the 
author used a Pi controller for correcting the error between the reference and measure current for estimating 
the stator resistance. This method has shown good performances and large stability, but using the two filters 
is able to produce a response delay in the system. In Figure 7 is performance at low speed and zero speed. 

The stability of the system response at low and zero speed has been checked by a rectangular 
command starting from zero to 300 rpm and with 50% of the rated load torque. by Figure 7(a) and 7(b) we 
have seen a very good stability and dynamic response compared to [8] who realized the same work but with 
a mechanical adaptation with a Pi regulator for the estimation of the stator resistance, which fell into 
a chattering problém at low speed with an average response time. On the other hand, we have noticed in 
the estimated and actual speed also the estimated parameters presents digital noises which are notable from 
the pulsations in the speed figures; this is mainly related to the incremental encoder and dSPACE card 
resolution [20]. Nevertheless, at zero reference, the sensorless drive shows a good stability. The torque 
presents good dynamics, but oscillations occur as it can be seen in the waveforms. 
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Figure 6. Performance of SM-DTC with speed estimation MRAS-PI with NN parameter estimation 
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Figure 7. Performance at low speed and zero speed 
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Figure 7. Performance at low speed and zero speed 


6. CONCLUSION 

In this article, we presented the implementation of a DTC control without a speed sensor using 
the rotor flux MRAS-based. Furthermore, we dealt with the resolution of the problems of parameters’ 
variations by an online estimation of the stator resistance and the inverse of rotor-time constant. We changed 
the conventional PI controller by two types of non-linear regulators; a sliding mode controller with a reaching 
law method in the speed control, and a neural network controller in the online parameters estimation. 
We validated our studies with the experimental results that have shown good dynamics in both high and low 
speeds, and even at zero speed reference. Our next work will focus on a more stable MRAS method, and on 
generalizing all the adaptation mechanisms by a nonlinear artificial neural network controller. 
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